The Approximate Capacity of 
the Gaussian A^-Relay Diamond Network 

Urs Niesen and Suhas N. Diggavi 



Abstract 

We consider the Gaussian "diamond" or parallel relay network, in which a source node transmits a message 
to a destination node with the help of N relays. Even for the symmetric setting, in which the channel gains to the 
relays are identical and the channel gains from the relays are identical, the capacity of this channel is unknown in 
general. The best known capacity approximation is up to an additive gap of order N bits and up to a multiplicative 
gap of order A^^, with both gaps independent of the channel gains. 

In this paper, we approximate the capacity of the symmetric Gaussian A^-relay diamond network up to an 
additive gap of 1.8 bits and up to a multiplicative gap of a factor 14. Both gaps are independent of the channel 
gains and, unlike the best previously known result, are also independent of the number of relays N in the network. 
Achievability is based on bursty amplify-and-forward, showing that this simple scheme is uniformly approximately 
optimal, both in the low-rate as well as in the high-rate regimes. The upper bound on capacity is based on a careful 
evaluation of the cut-set bound. We also present approximation results for the asymmetric Gaussian iV-relay diamond 
network. In particular, we show that bursty amplify-and-forward combined with optimal relay selection achieves 
a rate within a factor 0{\og*{N)) of capacity with pre-constant in the order notation independent of the channel 
gains. 

I. Introduction 

Cooperation is a key feature of wireless communication. A simple canonical channel model capturing 
this feature is the "diamond" or parallel relay network introduced by Schein and Gallager [[T]|, This 
network consists of a source node connected through a broadcast channel to N relays; the relays, in turn, 
are connected to the destination node through a multiple-access channel (see Fig. [T]). The objective is 
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Fig. 1. The A'^-relay diamond network. The source node u transmits a message to the destination node w via the N relays {v„}^^i. The 
two cuts indicated in the figure are the broadcast cut (separating the source u from the relays {i^n}) and the multiple-access cut (separating 
the relays {vn} from the destination w). 

to maximize the rate achievable between the source and the destination with the help of the N relays. 
Throughout this paper, we will be interested in the Gaussian version of this problem, in which both the 
broadcast and the multiple-access parts are subject to additive Gaussian noise. Moreover, for simplicity 
we will restrict attention in a significant part of the paper to the symmetric case, in which the channel 
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gains within the muhiple-access part and within the broadcast part of the network are identical (but are 
allowed to differ between the multiple-access and broadcast parts). However, we do show that some of 
the results for the symmetric setting can be extended to the asymmetric setting. 

For the Gaussian 2-relay diamond network, the rates achievable with decode-and-forward and with 
amplify-and-forward at the relays were analyzed in |l2l. It is shown there that these schemes achieve 
capacity in some regimes of signal-to-noise ratios (SNRs) of the broadcast and multiple-access parts of 
the diamond network. The asymptotic behavior of the A^-relay Gaussian diamond network was investigated 
in [[3]|. In certain regimes of SNRs of the broadcast and multiple-access parts of the network, it is shown 
that amplify-and-forward is capacity achieving in the limit as — > oo. New achievable schemes for 
the Gaussian diamond network with bandwidth mismatch (i.e., the source and the relays have different 
bandwidth) were introduced in dU and Perhaps surprisingly, these schemes lead to higher achievable 
rates than the ones obtained with amplify-and-forward and decode-and-forward even when the bandwidths 
at the source and the relays are identical. Half-duplex versions of the Gaussian diamond network, in which 
the relays cannot receive and transmit signals simultaneously, were considered in [6] and [7J. The capacity 
of a special class of 2-relay diamond networks is derived in [[8||. For networks in this class, one relay 
receives the signal sent at the source without noise, and the destination node is connected to the relays by 
two orthogonal bit pipes of fixed rate. To the best of our knowledge, this is the only non-trivial example 
for which the capacity of the diamond network is known for all values of SNR. For the general Gaussian 
iV-relay diamond network, the capacity is unknown. 

Given the difficulty of determining the capacity of communication networks in general and of the 
diamond network in particular, it is natural to ask if it can at least be approximated. For high rates, such 
an approximation should be additive in nature, i.e., we would like to determine capacity up to an additive 
gap. For low rates, such an approximation should be multiplicative, i.e., we would like to determine 
capacity up to a multiplicative gap. If a communication strategy can be shown to have both small additive 
as well as multiplicative gaps, then this strategy is provably close to optimal both in the high rate as well 
as low rate regimes. 

Additive approximations for channel capacity of communication networks were first derived in BU, 
where the capacity region of the two-user Gaussian interference channel is determined up to an additive 
gap of one bit. This was mainly enabled through a new outer bound for the interference channel. The 
approach of approximate capacity characterization was applied to general relay networks with single- 
source multicast in ifTOl . By introducing a new relaying strategy termed quantize-map-forward, capacity is 
derived up to an additive gap of Ibn bits, where n is the number of nodes in the network. This additive gap 
was improved through the use of vector quantization at the relays ifTTI . [|T2l . The sharpest known additive 
approximation gap is 1.26n bits for the complex Gaussian case (or 0.63n for the real case) [|T2ll . Since 
the A^-user diamond network is a special case of a relay network with a single source and destination 
and with n = + 2 nodes, these results yield an additive approximation up to a gap of 0.63A^ + 1.26 
bits for this network (assuming real channel gains). 

Multiplicative approximations were mostly analyzed for large wireless networks, for which the rate 
per source-destination pair is low. For a network with n nodes, the emphasis is on finding capacity 
approximations up to a small multiplicative factor in n. This approach was pioneered in [fT3l . Under a 
restricted model of communication, (essentially) the equal rate point of the capacity region of a wireless 
network with n randomly placed nodes was determined up to a constant multiplicative factor independent 
of n. Without the restrictive communication assumptions in |fT3l , the problem becomes considerably harder. 
Approximations for the equal rate point under a Gaussian model were derived in ffT4l up to a multiplicative 
factor of O(n^) for any e > 0. These approximation results were subsequently sharpened in ifTSl . lfT6ll to 
a factor n^^^/V^°s{n)) _ Uj^^jgj- gome conditions on the node placement, this factor can further be sharpened 
to 0(log(n)) (TT]. Multiplicative approximations for arbitrary relay networks with single-source multicast 
(as opposed to wireless networks with multiple unicast, i.e., multiple separate source-destination pairs) 
were derived in ifTOl . For a network with maximum degree d, the capacity is approximated to within a 
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factor of 2d{d + 1). As pointed out earlier, the Gaussian iV-relay diamond network is such a network 
with maximum degree d = N, and hence this result yields a multiplicative approximation up to a factor 

of 2iV(iV+ 1). 

To summarize, the capacity region of the general Gaussian A^-relay diamond network is not known. The 
best known additive approximation is up to a gap of 0.63A^+ 1.26 bits, and the best known multiplicative 
approximation is up to a factor of 2N{N + 1). In either case, the bounds degrade rather quickly as 
increases. It is hence of interest to find approximation guarantees that behave better as a function of the 
number of relays N in the network. Ideally, we would like the approximation guarantees to be uniform 
in in the network size. 

As a main result of this paper, we show that such a uniform approximation is indeed possible. More 
precisely, we find an additive approximation of the capacity of the symmetric Gaussian A^-relay diamond 
network of gap at most 1.8 bits for any SNR and number of relays A^. Moreover, we find a multiplicative 
approximation to the capacity up to at most a factor 14, again for any SNR and number of relays A^. This 
is a significant improvement over the previously best known additive approximation of 0.63A^ + 1.26 bits 
and multiplicative approximation of a factor 2N(N + 1), especially for large values of A^. In particular, 
as far as we know, this is the first such approximation result (both multiplicative as well as additive) that 
is independent of the number of network nodes for a nontrivial class of wireless networks. 

We further show that bursty amplify-and-forward (first introduced in [2, p. 76]) with properly chosen 
duty cycle is close to capacity achieving for the diamond network simultaneously in the sense of multi- 
plicative and additive approximation up to the aforementioned gaps. Hence, bursty amplify-and-forward 
with appropriately chosen duty cycle is a good communication scheme for the symmetric Gaussian A^-relay 
diamond network both at low and at high SNRs, and independently of the number of relays A^. 

Some of these results can be extended to the asymmetric setting. For general (i.e., not necessarily 
symmetric) Gaussian A^-relay diamond networks, we provide a factor 0(log'^(A^)) multiplicative approxi- 
mation of capacity, with pre-constant in the order notation independent of the channel gains. Achievability 
is based again on bursty amplify-and-forward, but this time a careful selection of relays is also necessary. 

The main technical contribution of this paper is the upper bound on capacity. The standard way to obtain 
upper bounds on the capacity of the diamond network is to evaluate two particular cuts in the wireless 
network, namely the one separating the source from the relays (called the broadcast cut in the following) 
and the one separating the relays from the destination (called the multiple-access cut in the following) 
as depicted in Fig. [T] This approach is taken, for example, in [|3l-[15]|. In fact, for symmetric Gaussian 
A^ -relay diamond networks, whenever the capacity is known, it coincides with the minimum of these two 
cuts. We show in this paper that, in order to obtain uniform additive or multiplicative approximations 
for the capacity of this network, considering just these two cuts is not sufficient. Instead we need to 
simultaneously optimize over all possible 2^ cuts separating the source from the destination. Without this 
careful outer bound evaluation, we believe that the uniform (in network size) approximation would not 
have been possible. 

The remainder of this paper is organized as follows. Section |ll] formally introduces the problem 
statement. Section [ni] presents the main results; the corresponding proofs are presented in Section HVl 
Section IVl contains concluding remarks. 

II. Problem Statement 

Consider the Gaussian A^-relay diamond network as depicted in Fig. |2l The source node u transmits a 
message to the destination node w with the help of A^ parallel relays {vi, . . . , vn}- The channel inputs at 
time t e IN at nodes u and f„ are denoted by X[t] and X„[t], respectively. The channel outputs at time 
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Fig. 2. The Gaussian A'^-relay diamond network. 



t e IN at nodes w and are denoted by Y[t] and Yn[t]. The channel inputs and outputs are related as 

Y4t]^y^X[t] + Z^[t], 

N 
n=l 

where (Z[t])t, ^ are independent and identically distributed Gaussian random variables with mean 

zero and variance one, independent of the channel inputs. The channel gains {gn)n=i {hn)n=i are 
assumed to be real positive numbers, constant as a function of time, and known throughout the network. 
A T-length block code for the diamond network is a collection of functions 

/: {1,...,M}^R^ 
/„:R^^R^, VnG{l,...,iV}, 
0: ^ {1,...,M}. 

The encoding function / maps the message W, assumed to be uniformly distributed over the set {1, ... , M}, 
to the channel inputs 

(^M)Li = f{w) 

at the source node u. The function /„ maps the channel outputs (l^nM)^i to the channel inputs 

at relay The decoding function maps the channel outputs at the destination node w into 

a reconstruction 

We say the code satisfies a unit average power constraint if 

if;E(X2[t])<l, 

t=i 

^Ee(X^M)<1, VnG{l,...,iV}. 

t=l 

'Note that the functions fn at the relays are not causal. This is to simplify notation; due to the layered nature of the network all results 
remain the same if causality is imposed. 
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The rate of the code is 

log(M)/T, 

and its average probability of error 

V{W ^ W). 

A rate R is achievable if there exists a sequence of T-length block codes with unit average power constraint 
and rate at least R such that the average probability of error approaches zero as T — t- oo. The capacity 
C(^N, {g„), (hn)) of the diamond network is the supremum of all achievable rates. 

A natural scheme for the diamond network is amplify-and-forward, in which each relay transmits a 
scaled version of the received signal. Formally, 

Xn[t] = anYn[t] = a„^/^X[t] + anZn[t], 

where the constant a„ is chosen to satisfy the power constraint at the relay. Denote by Ri (^N, (gn), (hn)) the 
rate achieved by amplify-and-forward with optimal choice of («n)^=i- We point out that the optimization 
over (an)n=i is necessary. While perhaps counterintuitive, it turns out that in the asymmetric setting the 
rate of amplify-and-forward is not always maximized when the relays use all available transmit power 
(see O for a discussion of this phenomenon in the two-relay case). 

If the SNR at the relays is low (i.e., gn <^ 1), it can be shown that simple amplify-and-forward is 
arbitrarily suboptimal. This is because the received signal power g^ at the relay f„ is much smaller than 
the noise power 1, and therefore the relay amplifies mostly noise. This effect can be mitigated by using 
bursty amplify-and-forward [;2|. For a constant 6 E (0, 1], called the duty cycle in the following, we 
communicate for a fraction S of time at average power 1/6 using the amplify-and-forward scheme and 
stay silent for the remaining time. This satisfies the overall average unit power constraint. The resulting 
achievable rate is denoted by Rs{N, (gn), (hn)) - This notation is consistent, i.e., for 5 = 1 the simple and 
bursty amplify-and-forward schemes coincide and achieve both rate Ri[N, (gn), (hn))- 

A special case of the general diamond network described so far is the symmetric setting, in which 
gi = g2 = ■ ■ ■ = gN = g and hi = h2 . . . = = h. With slight abuse of notation, we denote the 
capacity and rates achievable by bursty amplify-and-forward for the symmetric setting by C{N,g, h) and 
R6{N,g,h). 

Throughout this paper, we use bold font to denote vectors and matrices, log(-) and In(-) denote the 
logarithms to base 2 and e, respectively. All capacities and rates are expressed in bits per channel use. 



in. Main Results 

The main results of this paper are additive and multiplicative capacity approximations for the Gaussian 
diamond relay network. We start with a discussion of symmetric networks in Section IIII-AI General 
asymmetric networks are treated in Section IIII-BI 



A. Symmetric Diamond Networks 

The first result lower bounds the rate achievable over a symmetric diamond network by using bursty 
amplify-and-forward with optimized duty cycle 5. 

Theorem 1. For every symmetric diamond network with N > 2 relays and channel gains g,h > 0, there 
exists a duty cycle 6* G (0, 1] such that bursty amplify-and-forward achieves at least the rate 



( 1 

2 



Rs*{N,g,h) > { 



log (1 + lNmm{g,Nh}), 
|ln(4/3)log(l + iV^7), 
|log(l + iiV2^/i), 
iln(4/3)log(l + iVv^), 
[iln(4/3)log(l + iV2;^), 



ifmax{g,Nh} > 1 

ifma:x{g,Nh} < 1, g < h 

ifmax{g,Nh} < 1, g e {h,N^h), N^> 1 

ifmsix{g,Nh} <l, g E {h,N'^h), Ny/gh < 1 

ifma.x{g,Nh} < 1, g > N^h. 
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C{N,g,h) < { 



The proof of Theorem [T] is presented in Section IIV-AI Note that the optimal duty cycle 5* is allowed 
to depend on A^, g, and h. In the high-rate regime, i.e., the first and third cases in Theorem [H the duty 
cycle achieving the lower bound is 5* = 1, and hence the bursty amplify-and-forward scheme reduces to 
simple amplify-and-forward. On the other hand, in the low-rate regime, i.e., the second, fourth, and fifth 
cases in Theorem [U 5* < 1, and (genuine) bursty amplify-and-forward is used. 

Having established an achievable rate, the next theorem provides an upper bound on the capacity of 
the diamond network. 

Theorem 2. For every symmetric diamond network with N > 2 relays and channel gains g,h > 0, 
capacity is upper bounded by 

^ ^\og{l + Nmm{g,Nh}), ifmax{g,Nh} > 1 

l\og{l + Ng), ifmax{g,Nh}<l,g<h 

i log (1 + 2A^2^/i) +i, ifmax{g,Nh} < 1, g E {h,N^h), Ny/^> 1 

log(l + 2iVv^), ifma.x{g,Nh} < I, g e {h,N^h), Ny/^< 1 

^\\og{l + N^h), ifmax{g,Nh} < 1, g > N^h. 

The proof of Theorem [2] is presented in Section IIV-BI As a corollary to Theorems [T] and [21 we obtain 
that bursty amplify-and-forward is close to optimal, in the sense that it achieves capacity both up to a 
constant additive gap as well as a constant multiplicative gap, where both constants are independent of the 
number of relays and the channel gains g and h. This shows that optimized bursty amplify-and-forward 
is a good communication scheme for the symmetric diamond network both at low rates (due to the small 
multiplicative gap) as well as at high rates (due to the small additive gap). 

Corollary 3. For every symmetric diamond network with N > 2 relays and channel gains g,h > 0, there 
exists a duty cycle 6* G (0, 1] such that 

C{N,g,h) - R5*{N,g,h) < 1 + | log(3) < 1.8 bits, 

and 

C{N,g,h) 



Rs*{N,g,h) - ln(4/3) 

The proof of Corollary [3] is presented in Section IIV-CI We point out that choosing the duty cycle S* as 
a function of N, g, and h, is not necessary to obtain the additive approximation result in Corollary [3l In 
fact, using only simple amplify-and-forward achieves the same additive approximation guarantee, i.e., 

C{N,g,h) - Ri{N,g,h) < 1.8 bits 

for all N > 2, g,h > 0. However, the same is not true if we are also interested in multiplicative 
approximation guarantees (at least in the low-rate regime). To achieve a constant additive approximation 
as well as constant multiplicative approximation, the duty cycle 6* is required to vary as a function of N, 
g, and h, and therefore bursty amplify-and-forward is required. 

From Theorems [T] and [2l the capacity of the symmetric diamond network has three distinct regimes, 
depending on whether g < h, h < g < N'^h, or g > N'^h. In the first regime {g < h), the channel gain to 
the relays is weak compared to the channel gain to the destination, and the achievable rate is constrained 
by the broadcast part of the diamond network. The capacity in this regime is given approximately by 

C{N,g,h)^l\og{l + Ng), 

where the approximation is in the sense of Corollary [3l namely up to a multiplicative gap of factor 14 in 
the low-rate regime (g <^ N^^) and up to an additive gap of 1.8 bits in the high-rate regime (g ^ A^~^). 
This is the capacity of a single-input multiple-output channel with unit power constraint, one transmit 
antenna, N receive antennas, and channel gain ^/g between each of them. Thus, the broadcast cut in 
Fig. [T] in Section U is approximately tight in this regime. 
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In the third regime {g > N'^h), the channel gain to the relays is strong compared to the channel gain 
to the destination, and the achievable rate is now constrained by the multiple-access part of the channel. 
The capacity in the third regime is given approximately by 

C{N,g,h)^'^\og{l + N^h). 

This is the capacity of a multiple-input single-output channel with unit per-antenna power constraint, 
transmit antennas, one receive antenna, and channel gain \/h between each of them. Thus, the multiple- 
access cut in Fig. \T\ is approximately tight in this regime. Observe that to achieve this rate the signals 
sent by the relays must be highly correlated and add up coherently at the destination. 

The most interesting regime is the second one (h < g < N'^h). If ma,x{g, Nh} > 1, then the capacity 
is given approximately by 

C{N, g,h)^\ log {l + N min{(7, Nh}) , 

and again either the broadcast cut or the multiple-access cut are tight. If ma,x{g, Nh} < 1 the situation 
is more complicated. If N^/gh > 1, then the capacity of the diamond network is approximately 

C{N,g,h)^l\og{l + N'gh), 

and, if N^/gh < 1, 

C{N,g,h) ^\\og {l + N^). 

In both cases, the capacity depends on the product of g and h, and not merely on the minimum of g 
and Nh. Hence, neither the broadcast cut nor the multiple-access cut are tight in this case. In fact, these 
bounds can be arbitrarily bad, both in terms of additive gap as well as multiplicative gap, as the next two 
examples illustrate. 

For the additive gap, consider g = N^^/^ and h = N~^^^. Then max{g, Nh} = N^'^^^ < 1, g = 
N^/^h G {h, N^h), and Ny/^ = N^/^ > 1, so that 

C{N,g,h)^'^\og{l + N^gh) 
= |log(l + iVi/^). 

On the other hand, the minimum of the broadcast and multiple-access cuts yields 

I log {l + Nmm{g,Nh}) = i log(l + iV=^/«), 

resulting in an additive gap of order 6(log(7V)) bits, which is unbounded as the number of relays N oo. 

For the multiplicative gap, consider g = A^~^ and h = N~^. Then maxj^f, Nh} = N^"^ < 1, g = Nh E 
{h, N^h), and N^ = A^-^/^ < 1, so that 

C{N,g,h)^l\og{l + N,/^) 
= ilog(l + iV-3/2) 

^ |log(e)iV-=^/2_ 

On the other hand, the minimum of the broadcast and multiple-access cuts yields 

I log (1 + Nmm{g, Nh}) = | log(l + N'^) 

^|log(e)iV-\ 

resulting in a multiplicative gap of order 6(v^), which is again unbounded as the number of relays 
oo. 

In the second regime, we thus need to take cuts other than the broadcast and multiple-access ones into 
account. The need for this can be understood as follows. Consider a general cut separating the source node 
u from the destination node w in the diamond network as shown in Fig. [3l Formally, let 5 C {1, . . . , A^}, 
and consider the cut from u U {f„}„e5 to u; U {t'„}„e5<:. Assume the signals {Xn)n=i sent from the relays 
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to the destination are highly correlated. This results in the signal summing up coherently at the receiver, 
increasing the rate across the cut. At the same time, if the signals sent from the relays are highly correlated, 
then the signals (X„)„e5<= available at the relays on the other side of the cut can be used to estimate the 
signal received at the destination node. This decreases the rate across the cut. Thus, for general cuts, 
there is a tradeoff between the gain from coherent reception and the loss from prediction that come with 
increased signal correlation. This tradeoff is absent if we only consider the broadcast and multiple-access 
cuts. It is precisely this tradeoff that determines the behavior of the capacity of the diamond network in 
the second regime. 




Fig. 3. A general cut in the diamond network. Here 5 C {1, . . . , A''}, and the cut separates u U {unjnes from w U {^nlnes^- 

We point out that a (partial) decode-and-forward strategy is not sufficient to provide a uniform capacity 
approximation as in Corollary [3l Indeed, due to symmetry, all relays would be able to decode the source 
in any such strategy, which implies that decode-and-forward and partial decode-and-forward coincide in 
this case. The rate achievable with decode-and-forward is given by 

ilog (1 + mm{g,N^h}). 

Comparing this with Corollary [3l we see that (partial) decode-and-forward has an additive gap of at least 
f2(log(A^)) bits and a multiplicative gap of at least a factor Q(N) to capacity. Similarly, compress-and- 
forward does not achieve constant (in the network size A^) additive or multiplicative gaps to capacity, 
since it does not capture the gain from coherent signal addition at the destination. Finally, as was pointed 
out earlier, the traditional amplify-and-forward strategy does not yield a constant factor approximation of 
capacity. In fact, in can be shown that simple amplify-and-forward results in unbounded multiplicative 
gap even for = 2. Therefore the bursty amplify-and-forward scheme introduced in ^ and advocated 
in this work has the nice property of being uniformly approximately optimal in both the additive sense 
and the multiplicative sense, as well as being a simple modification of the traditional amplify-and-forward 
scheme. 

B. Asymmetric Diamond Networks 

In the last section, we have considered symmetric diamond networks, in which the channel gain from 
the source u to the relay Vn is ^ and the channel gain from Vn to the destination w is \/h for all n. 
In this section, we show how some of the results can be extended to asymmetric diamond networks, in 
which the channel gains {gn)n=i ^^^^ {hn)n=i are allowed to take arbitrary values. 

For this asymmetric setting, it was shown in IfTOl that (bursty) amplify-and-forward does not achieve a 
constant (as a function of Qn and /i„) additive-gap approximation even when N = 2. However, we show 
here that bursty amplify-and-forward is approximately optimal in the sense of multiplicative approximation 
for any A^, {gn)n=i^ and {hn)^^^. More precisely, we show that the rate achieved by bursty amplify-and- 
forward combined with optimal relay selection is at most a factor 0(log^(A^)) from capacity uniformly 
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in {gn)n=i ^iid {hn)n=i- While not constant in A^, compared to the best previously known multiplicative 
approximation of a factor 2N{N + 1), this is a significant improvement. Hence, at least in the low- 
rate regime, bursty amplify-and-forward is also a good communication scheme for asymmetric diamond 
networks. 

Theorem 4. There exists a universal constant K < oo such that for every diamond network with N >2 
relays and channel gains {gn)n=i^ (^n)^=i > 0, 

V 'Vy'^^'V njj <K\og\N). 



The proof of Theorem |4] is presented in Section ITV-Di At a high level, achievability is proved as follows. 
Group the relays into classes such that all relays in the same class have approximately the same channel 
gains. Choose one such class, and set the constants a„ = for all relays not in this class (effectively 
disabling those relays). This relay-selection step reduces the original asymmetric network to a (almost) 
symmetric subnetwork. Theorem \T\ can now be applied to this symmetric subnetwork to obtain a lower 
bound on the rate achievable with bursty amplify-and-forward. By maximizing over all possible classes, 
and hence all possible symmetric subnetworks, we get the largest rate achievable in this manner. The 
corresponding upper bound shows that this approach of relay selection combined with bursty amplify- 
and-forward is approximately optimal. 

IV. Proofs 

Sections IIV-A[ IIV-B[ and IIV-CI contain the proofs of Theorem \T\ (achievability). Theorem [2] (upper 
bound), and Corollary |3] (approximation) for symmetric diamond networks. Section IIV-DI contains the 
proof of Theorem |4] for general asymmetric diamond networks. 

A. Proof of Theorem [7] (Achievability for Symmetric Networks) 

We start with a lemma computing the rate achievable by amplify-and-forward. 

Lemma 5. For every symmetric diamond network with N > 2 relays and channel gains g,h > 0, 
amplify-and-forward achieves 

N^gh 



i?i(iV,^7,/i) = ilog(l + - 



+ g + Nh 

Proof: Recall that with amplify-and-forward relay Vn transmits 

Xn[t] = aYn[t] = a^X[t] + aZ^t] 

at time t, with constant a > chosen to satisfy the average unit power constraint. The received signal at 
the destination node w is 

N 



Y[t] = aNy^X[t] + ay/hJ2^n[t] + Z[t]. 



n=l 

Observe that this describes a memoryless point-to-point channel with channel gain aN-/gh and additive 
Gaussian noise of variance 1 + a^Nh. Ri{N, g, h) is the capacity of this point-to-point channel, optimized 
over the value of a. 

For any value of a > 0, the optimal distribution of the input X for this point-to-point channel is 
Gaussian with zero mean and variance one. The signal sent by the relays has power 

nxl) = a\l + g), 
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and hence for 

a'e [0,1/(1 + (7 

the average unit power constraints at the relays are satisfied. This implies that amplify-and-forward achieves 
a rate of 



Ri {N, g, h) = max i log (l + " ) 

N^gh 



i log 1 



2 



1 + g + Nh 



The next lemma describes the rate achievable with bursty amplify-and-forward. 

Lemma 6. For every symmetric diamond network with N > 2 relays and channel gains g,h > 0, bursty 
amplify-and-forward with duty cycle 6 G (0, 1] achieves 

Proof: During the 5 fraction of time we communicate, we are dealing with an equivalent channel 
with gains VgJS, Vh/S and with unit power constraints. The result now follows from Lemma |5]by taking 
into account that we only communicate a fraction S of time. ■ 
Note that Lemmas [5] and [6] coincide for 5 = 1, as expected. We now proceed to the proof of Theorem [U 
To simplify notation, set 

Rs^Rs{N,g,h) 

for 6 e (0, 1]. 

We consider the cases raax{g, Nh} > 1 and maxj^f, Nh} < 1 separately. Assume first maxj^f, Nh} > 
1. Here we set 6=1, i.e., we use simple amplify-and-forward. By Lemma [5] 



R 



N^gh 

' = -^'''^^'^^+g + Nh 



>g (l + - 



minj^f, Nh} maxjgf, Nh} 



minj^f, Nh} + maxj^f, Nh} 
N min{g, Nh} max{g, Nh} 



^ ' ' 3 max{g, Nh} 

= l\og{l + lNmm{g,Nh}), 

where we have used that 1 < max{g, Nh} to obtain the inequality. 

Assume in the following that max{g, Nh} < 1. We consider the cases g < h, g E {h,N'^h), and 
g > N'^h separately. Consider first g < h. Bursty amplify-and-forward with duty cycle 5 = Ng < Nh < 1 
achieves by Lemma |6] 

N'gh/{N'g'^ 



Rs = \Ng log (l + y^iy^^^^) \ 

2 y l + g/{Ng) + Nh/{Ng)) 



\Ng\og (l + 



h 



g + g/N + h 



W . , / h 



(a) / 

> lNg\og[l 



h + h/N + h 

> 'jNg \og{4/3) 

> iln(4/3) log(l + A^^?), 
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where in (a) we used g < h. 

Consider then g G {h,N'^h). If Ny/gh > 1, then simple amplify- and-forward achieves by Lemma [5] 

>ilog(l + |iVV), 

where we have used that 1 + g + Nh < 3, which follows from ma,x{g, Nh} < 1. 

Still assuming g e {h,N'^h), if Ny/gh < lH then bursty amplify-and-forward with duty cycle 5 -- 
Ny^gh < 1 achieves by Lemma [6] 



Rs = \N^\og (l + - 
iiVv^log (l + - 



N^ghl{N^gh) 



g/{N^) + Nh/{N^) 



1 



2'-^^--^y ' 1 + v^/(iVv^) + v^/v^ 

= |iVv^log(4/3) 
>iln(4/3)log(l + iVv/^), 

where in [h] we have used that g < N'^h and g > h. 

Consider finally g > N'^h. Bursty amplify-and-forward with duty cycle 5 = N'^h < g < I achieves by 
Lemma [6] 

N^gh/{N^h^] 



'^'='^'''^^'4'^l + g/iN^h) + Nh/iN^h) 
^iiV^Mog(l + - « 



g/{N^h) + 1/N 



(c) 



> lN'h\og{4/3) 



> |ln(4/3) \og{l + N^h), 
where in (c) we have used that < 1 < g/{N'^h). 



B. Proof of Theorem |2] (Upper Bound for Symmetric Networks) 

In this section, we derive an upper bound on the capacity of the Gaussian diamond network. The 
standard way to find such bounds is to start with the cut-set bound and then to simplify it further to 
obtain a closed-form expression. The derivation here starts with the cut-set bound as well, but differs in 
several key aspects from the standard approach, which we now highlight. 

Let 

[iV]4{l,2,...,iV}, 

and for a subset S (Z [N], define 

S'' 4 [N] \ S. 

By the cut-set bound [[HI Theorem 14.10.1], 

C{N,g,h) < sup min l{X,Xs;Y,Ys^ \ Xs^), 

^Note that g G {h, N^h) and N\^i < 1 imply max{g, Nh} < 1. 



12 



where the maximization is over random variables X.Xy^] satisfying the power constraints E(X^) < 1, 
E(X2) < 1, and where Xg = for any subset S C [N] (see Fig. [3] in Section HlFAl) . A short 

calculation (done in Q below) reveals that 

sup min l{X,Xs;Y,Ysc \ X^c) < sup min (l{X;Ysc) + l{Xs;Y \ XsM. (1) 

In the right-hand side of O, the first mutual information corresponds to the rate between the source 
nodes and the relays, and the second mutual information corresponds to the rate between the relays and 
the destination node. 

One approach is to simplify this expression further through a sequence of two steps. The first step is 
to upper bound 



l{Xs; Y I Xsc) =n{Y\ Xs^) - n[z) 

<n{VhEnesXn + z)-niz) 



where, in order to avoid confusion with the channel gain h, we denote the differential entropy by the 
non-standard symbol "H. This first step thus removes the conditioning on the signals Xgc available at the 
destination side of the cut. The second step is to interchange the order of maximization and minimization. 
This yields 

C{N,g,h)< mm sup (/(X; F^O + v^Enes^n + ^ 

(i log(l + {N- n)g) + i log(l + n'h)) . (2) 



= min 

ne{0,...,Af} 

This can be further upper bounded by considering only = or n = iV, resulting in the minimum of 
the broadcast and multiple-access cut 

C(iV, h) < min {i log(l + Ng),\ log(l + N^h)]. (3) 

Neither of the upper bounds ^ and ([3]) are tight enough to obtain a constant gap approximation of the 
capacity (this can be seen from the two examples presented after Corollary [3]). 

In this paper, we also start the derivation of the upper bound from the cut-set bound ([1]), but we avoid 
taking the two simplifying steps mentioned in the last paragraph. Instead, we first show, using the symmetry 
in the problem, that the correlation between any two signals X„ and X^ with n ^ h can be assumed to 
be equal without loss of optimality. Using the resulting simple form of the covariance matrix allows us 
then to evaluate the term /(X^; Y\Xsc) directly. This enables us to keep the conditioning on X^c, which 
yields a significantly tighter upper bound on capacity. The resulting upper bound is summarized in the 
following lemma. 

Lemma 7. For every symmetric diamond network with N > 2 relays and channel gains g,h > 0, capacity 
is upper bounded as 

C{N,g,h) 

< sup min ( Uog{l + {N - n)g) + Uogfl + n(l + {n - l)p - — )h)\ 

The variable p appearing in the lemma can be interpreted as the correlation between the random variables 
X[Ar] as mentioned in the preceding discussion. Note that it is not clear a priori that this correlation p can 
be restricted to be nonnegative. This restriction is part of the assertion of the lemma. We also point out 
that it is important that p = 1 is excluded from the supremum in Lemma |7l the result is not true without 
this restriction. 
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It will be convenient in the following to work with a weaker version of Lemma U\ Note that, for 

pe [0,1), 

1 + (n - l)p )— < 1 + np ^ 



{N -n-l)p ~ l + {N-n)p 

N \f^ + {N-n)p 



< 



N-nJ\ l + {N-n)p 
N 



N -n 



Hence 



C(iV, g,h)< min f | log(l + {N - n)g) + \ log (l + ^^h) \ . (4) 

The upper bound (Hj) derived from Lemma |7] can be compared to the simpler bound If n = KN 
for some constant K E (0, 1), then the factor multiplying the channel gain /i in (|2l) is of order 6(A^^). On 
the other hand, the same factor in Q is of order Q(N). Thus, the bound dH) can be considerably tighter 
than the simpler bound 

Proof of Lemma ^ By the cut-set bound [[TSl Theorem 14.10.1], 

C ^ CiN, g, h) < sup min J(X, Xs; Y, Ys^ \ Xs^) , (5) 

where, as before, the maximization is over random variables X, X[7v] satisfying the power constraints 
E(X^) < 1, E(X^) < 1. We evaluate dD in two steps. First, we argue that the maximization over X, X[n] 
can be restricted to jointly Gaussian random variables such that each E(X^) = 1 and 'E(XnXn) = p for 
n ^ h and some p E [— 1/(A^ — 1), 1]. This simplifies the maximization to be over just the parameter p 
instead of A^-dimensional distributions. Second, using the resulting simple form of the input distributions, 
we analytically evaluate the mutual information in ^ to obtain the stated bound. 

We start by simplifying the mutual information in ^ for a fixed cut S C [N] . We have 

l[X,Xs]Y,Ysc I Xgc) 

= n (F, 1^5= I Xsc) - n {Y, Ysc I X, X[N]) 

= n{Ysc I Xsc)+n{Y I Ys^,Xsc) -H{Ysc \ X,X^m])-'H{Y \ Ys^,X,X^n^) 
<H{Ys.)+H{Y I Xsc) -H{Ys. \ X)-'H{Y \ X[^]) 

= l{X;Ys^)+l{Xs;Y\Xsc), (6) 
where we have used that 

n{Ysc I = n{Zsc) = u{Ysc \ x), 

and that 

U{Y I Ys^,X,X^M]) =H{Z)=H{Y \ X[n]). 
Combining © and ^ yields 

C < sup min (/(X; Ys^) + l{Xs; Y \ Xsc)). (7) 

For the first term in dV]), 

l{X;Ysc) <l\og{l + \S''\g), (8) 
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since the channel from X to Igc is a Gaussian single-input multiple-output channel with channel gains 
y^. For the second term in 

l{Xs; Y I Xsc) ='H{Y\ Xs^) - ^/(y | X[^]) 

< n{VhT.n^s{Xn - Pn{Xs^)) + - H{Z), (9) 

for any choice of functions f3n{Xsc) for n E S.ln particular, let f3n{Xsc) be the minimum mean-square 
error estimator for X„ based on X^c 

Let X[Ar] have covariance matrix Q. Then, by lfT9l Theorem 1.2.11], (X„ — /3„(X5c))„g5 has covariance 
matrix 

Qs\S'' = Qs,s — Qs,S''Qs'^,S''Qs'',s, (10) 

where, for any subsets 5*1, 5*2 C [N], Qsi,S2 is the submatrix of Q induced by the rows 5*1 and columns 5*2, 
and where Q^c 5= is the Moore-Penrose generalized inverse of the matrix Qs'=,S'=- The matrix Qs\s'' is called 
the generalized Schur complement of Qs'^^S'^ in Q. Note that if Qs'^^S'^ is invertible, then Q^c gc = Q^c^^c 
and the generalized Schur complement reduces to the standard Schur complement. 

Before proceeding, we need to introduce some notation. Denote by the ax a identity matrix, and by 
la b the a X b matrix of ones. To simplify notation, we will write 1 for the column vector 1^ 1, whenever 
the dimension is clear from the context. With these definitions, 

^(v^Enes(^n - f3n{Xsc)) + z) - < | log (l + hl^QslS^l). (11) 

Substituting ([8]), ©, and ^ into © yields 

C< sup min U\ogil + \S'\g) + l\og{l + hl^Qs\s^l)), 

Q>0: SclN] \ / 
'?n,n<lVne[7V] 

where Q > denotes that Q is a positive semi-definite matrix. We have thus simplified the maximization 
over input distributions to a maximization over covariance matrices. The next step is to show that the 
covariance matrix Q can be restricted without loss of optimality to have the form 



plAr,7V + (1 - p)/ 



and hence the maximization over covariance matrices can be further simplified to a maximization over 
just the scalar correlation parameter pJl 
For convenience of notation, define 

Vs(Q) = |log(l + W'^Qs|5=l) 

and 

= ( 5log(l + 1^1^?) + MQ) 

so that 



SclN] ^ 2 



C < sup i^{Q). (12) 

Q>0: 

q„,n<iyn£[N] 

'Upon completion of this work, we realized that a somewhat similar argument as in this step was used in |20t Section III] for the Gaussian 
multiple-access channel with feedback. 
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Consider a covariance matrix Q > 0, and let P be any permutation matrix on [N]. Note that P^QP > 0. 
Moreover, by symmetry]^ 

ijiQ) = min (ilog(l + \S^\g) + i,s{Q)) 
sg[n] 

= min (ilog(l + \S'\g) + tPs{P^QP)) 

sg[n] 

= i;{P^QP), 

and thus tp{-) is invariant under permutation. 

Now, the generalized Schur complement is matrix-concave over the set of positive semi-definite matrices 
[[2T1 Theorem 3.1] (see also [|22l p. 469] for the corresponding result for positive definite matrices). More 
precisely, if Q = XiQ^ + X2Q'^ with Ai G [0, 1], A2 = 1 - Ai, then 

i.e., 

Qs\S'= — i^lQsiS'' + -^zQlisc) 
is a positive semi-definite matrix. Therefore, 

implying that 

1 log (1 + hl^Qsis^l) > I log (1 + X.hl^Ql^s^l + \2hl^Q%s^l) 

> Aii log (1 + /il^Q^|5cl) + X2I log (1 + /il^Qli^cl). 
Thus ipsiQ) is concave in Q. Finally, 

mm (Ihgil + \S^\g) + MQ)) 

> min (Ai(ilog(l + \S'\g)+MQ')) + A2(|log(l + \S'\g) + MQ' 

> Ai min (ilog(l + \S'^\g) + MQ')) + A2 min (|log(l + \S'^\g) + MQ')), 

SC[N] ' SC\N\ ' 

and hence ^/'(Q) is also concave in Q. 

Fix £ > 0, and assume that Q* achieves er-optimality, i.e., Q* > 0, „ < 1 for all n G [A^], and 

^(Q'') > sup 

Q>0; 

<lVne[A'] 

Set 

P 

where the sum is over all A^! permutation matrices on [A^]. 

Note that Q is positive semi-definite and satisfies g„ „ < 1 for all n G [A^]. Moreover, using the 
concavity and invariance under permutation of we obtain 



p 

= mi. 

''Note that the minimization over 5 C \N\ is crucial for this fact to hold. Indeed, ipsiQ) 7^ ^s{P'^QP) in general. 
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and hence Q is also an e-optimal covariance matrix. Note that this Q has the form 

for K < 1 — p], and thus we can restrict the maximization of 'tp(Q) to matrices of this form. Since the 
generalized Schur complement is monotonically increasing over the set of positive semi-definite matrices 
[|2T1 Theorem 3.1], we can further restrict the value of k to be 1 — p. Denote the resulting matrix by Q'^, 
i.e.. 





/I 


P 


P ■ 






P 


1 


P ■ 


■ P 




P 


P 


1 . 


■ P 






P 


P ■ 


■ 1/ 



Note that is positive semi-definite only if p G [— 1/(A^ — 1), 1] (since otherwise the eigenvalue 
corresponding to the eigenvector 1 is negative). 

The upper bound on capacity in (fT2l) can thus be simplified to 



(13) 



C < sup min (llog(l + \S''\g) + i log (l + /il^Q^i^.l)) 



where [0] is understood as the empty set and [0]^ = [N]. Observe that the minimization in (fT3l) is over 
integers n G {0, . . . , A^} as opposed to subsets S C [A^] due to the symmetry in Q^. Note furthermore 
that instead of maximizing over arbitrary input distributions, we only have to maximize over the single 
real number p. 

We now compute the expression in parentheses in (fT3l) analytically. To this end, we need to compute 
Q[n]\[n]" which, by (flOl) . involves the computation of the generalized inverse • will first 

consider the case when Qj'^j^ ^^-^^ is invertible, and then consider the remaining cases in which Q^^j^ j^^^ is 
not invertible. If n G {1, . . . , — 1} and p G [— 1/(A^ — 1), 1), then Qj^^^c [„]c is invertible, and after some 
algebra, we obtain 



^[n]\[n]- V ^ ^' l + (N -n-l)pJ 

We now consider the remaining cases, in which Qj'^jc [„]c is not invertible. If p = 1 and n G {1, 
1}, then 

1 Q\n]\\nV''^ = 0- 



If n = 0, then 
and if n = A^, then 



mm 



I'^n^ 1 = 
.l = l^Qn = Nil + iN -l)p), 



(14) 
N- 
(15) 

(16) 
(17) 



both for any p G [-1/(A^ - 1), 1]. 
Denote by 

r]{p,n) = n[l + {n - l)p 



n{N — n)p'^ 
1 + {N -n-l)p. 

the right-hand side of GH). Note that r]{p,n) is well defined for all p G [-1/(A^ - 1), 1], n G {0, 
except for r](l,N) and 'i]{—l/(N — 1),0) (for which the expression involves dividing zero by zero). 



N} 
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Moreover, from (fT4l)-(fT71) we see that whenever r](p,n) is well defined, it is equal to l^Qj'^ju^j,,!- For 
the two cases in which r]{p, n) is not well defined, we have from (fT5l)-(fT7l) that for any n E [N], 



limr/(p,n) = 1 Q[„]|[„]cl, (18a) 

.,-1)^^-1)^^^'") = ^"^Mlff"^!' (18b) 

and in particular this holds for = and n = N. Thus we can write l-^Qj^^ju^j,,! compactly as a 
function of r]{p, n) for any n G [A^] and do not need to consider the boundary cases n = {i, n = N , and 
p = -1/{N - 1), p = 1 separately. Substituting ([H]) and into we obtain 

C< sup min (llogil + {N - n)g) + ^\og (l + n(l + {n - l)p - ^ 

(19) 

Observe that the supremum in (fT9l) is only over p G (— 1/(A^ — 1), 1) as opposed to p G [— 1/(A^ — 1), 1] 
as in (fT3]) . 

We finally argue that the supremum can be restricted to be over values p G [0, 1). Consider the derivative 
with respect to p of the multiplier of the h term in (fT9l ), 

—n 1 + (n - l)p ) ■^—r- = n (n - 1) - n(A^ - n)p ^ ^ . 

dp V ^ l + (iV-n-l)py ^ ^ ^^(l + (Ar_^_l)p)V 

If p G (— 1/(A^ — 1),0), then this derivative is non-negative, and thus the multiplier of h in (fT9l) is non- 
decreasing in that range of p. Since this is true simultaneously for all n G {0, . . . , A^}, we can restrict the 
supremum to be over the range p G [0, 1). This proves the lemma. ■ 
We now proceed to the proof of Theorem [21 As before, we denote the capacity of the diamond network 

by 

We again consider the cases max{(7, NK\ > 1 and maxj^f, Nh} < 1 separately. Assume first that 
max{g, Nh} > 1. Capacity is upper bounded by the minimum of the simple broadcast and multiple- 
access cuts 

C < min {i log(l + AT^?), i log(l + iV2^) } 

= ilog (1 + Armin{^,A^/i}). (20) 

Observe that ( |20l) is valid regardless of the value of maxj^f, Nh}. 

Assume in the following that raa,x{g, Nh} < 1. As before, we treat the cases g < h, g E {h,N'^h), 
and g > N'^h separately. Consider first g < h. Using the upper bound in (|20l ). we obtain 

C <l\og{l + Ng). 

Consider then g E {h, N'^h). If Ny/gh > 1, then the simplified form © of Lemma |7] with N — n = 
\N'^h'\ G {0, . . . , A^} (since Nh < 1) yields 

C < \ log(l + (A^ - n)g) + \ log (l + ^ 

f N^ \ 
= I log(l + \N'h]g) + i log i^l + J^hj 

< l\og{l+ g + N^gh) + l 
<ilog(l + 2ArV) + i 
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where we have used that g < N'^gh since 

N^h > V^\^Wh = iVv^ > 1. 
Still assuming g E {h, N'^h), if Ny/gh < 1, then the simplified form (HI) of Lemma |7] with N — n 



\NyJh/gA^ G {0, . . . , A^} (since g > h and hence \/h/g < 1) shows that 

C < i log(l + {N- n)g) + \ log + -/i 

= \\og{l+\N^^g) + \\og' ' 



< ilog (1 + ^7 + N^) + ilog (1 + N^) 
<log(l + 2iVv/^), 

where we have used that 

g < VgVWh = NVgh. 
Finally, consider g > N'^h. The upper bound (|20l) yields 

C< ilog(l + Ar2/i), 

concluding the proof. 

C. Proof of Corollary \3\ ( Capacity Approximation for Symmetric Networks) 
The corollary follows directly from Theorems [T] and [2] using the inequalities 

>alog(l + x), for a G [0, > 



log(l + ax 

and 



<alog(l + a;), fora>l,a;>0, 



log(l + ax) 



> log(a) + log(l for a G [0, > 
< log(a) + log(l + x), for a > 1, X > 0. 



D. Proof of Theorem |?] ( Capacity Approximation for Asymmetric Networks) 

The idea of the proof is as follows. Group the relays into classes such that all relays in the same class 
have approximately the same channel gains. We argue that the number of classes needed is on the order 
9(log^(A^)). Choose one such class, and set the constants a„ = for all relays not in this class. This 
effectively reduces the network to a (almost) symmetric one, which we have already analyzed in the earlier 
parts of this paper. By maximizing over which class to choose, we get the largest rate achievable in this 
manner. This yields a lower bound on Rs{N, (gn), {hn})- We then argue that this approach is close to 
optimal, by showing that capacity C(^N, (gn), {hn)) is upper bounded by 6(log'^(iV)) times the maximum 
of the capacities of these classes. 

Recall the notation 

[iV]^{l,...,iV} 

and, for S C [A^], 

Furthermore, in this section, we will use 

9s — (5'n)neS! 

hs — {hn)nes 
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for 5 C [iV], and 

ahs = {ahn)nes 

for scalar a G R. 

We want to partition [N] into subsets such that for n and h in the same subset the relays Vn and Vn have 
approximately the same channel gains. Moreover, we want the number of required subsets to be small. 
This is not directly possible if the channel gains are very different. For example, consider g-n = hn = 2"; 
note, however, that in this case most of the relays are very weak compared to the strongest one, and could 
hence be disregarded without too much loss in rate. We formalize this idea by allowing some "overload" 
subsets (in the language of quantization theory) in the partition of [A^], which correspond to relays that 
may have very different channel gains, but that are all too weak to have much impact on achievable rates. 

Define 

g* = max mm{gn, N'^hn}, 

nG[N] 

h* = max mini /in, fif„|. 

ne[N] 

The quantities g* and h* are essentially the largest channel gains, accounting for situations in which one 
of the channel gains hn clearly dominates the other one. If we let n be such that h* = min{/i„, 
then 

g* > mm{gn, N'^K} > ram{gn, h^} = h*. (21) 

Similarly, if n is such that g* = min{(7„, N'^hn}, then 

h* > mm{K, gn} > iV"' mm{N^K, g^} = N^^g*. (22) 

Thus, g* and h* can not be too different. 

We are now ready to introduce the partition of [N] mentioned above. We start with the "overload" 
subsets. Define the sets 

T'^{ne[N]:g^<N~'g*}, 
T^A [nE [N] \T' ■.h^< N-^h*}, 

i.e., and correspond to those relays that have channel gains that are very weak compared to the 
strongest one in the network. Set 

L4 L31og(iV)J. 

For £ e {0, ... , L}, define 

T/ ^{ne [N] \ {T' U T') : g^ G i2~'~'g\ 2"V], /i„ > }, 

r| ^ {ne [N] \ {T' U T' Tj) : g,, > N'hn, K E {2-'-^h\ 2-'h^] ] , 

i.e., {T/} and {T^} quantize those channel gains for which one of hn dominates the other one. Finally, 
define for k,£ E {0, . . . , L}, 

Sk,i = {ne [N] \ {T'UT' Ui{T}uT?)) : g^ G {2~'-'g\2-'g% G (2-^-l/^^ 2^^]}. 

The subsets {Sk/} quantize the remaining channel gains. The number of sets T^, T^, {T/}, {T/}, {Sk,e} 
is equal to 

L ^ (L + 1)2 + 2(L + 1) + 2 = Q{\og\N)). 

We argue that T^, T^, {T/}, {T/}, {Sk,e} partition [N]. The sets are clearly disjoint, so we only need 
to show that their union covers [A^]. If either g^ < N~^g* or /i„ < N^^h* then n G U T^. Assume 
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in the following discussion that Qn > N ^g* and hn > N ^h*. If < g* and /i„ < h*, then n is an 
element of {T/}, {T/}, or {5^,^. If g,, > g\ then 

/in < N^hn <g* < gn, 

SO that gn > N'^hn and hn = min{/i„, (7^} < /i*. This implies that n E UiT^. If /i„ > /i*, then 

g-a <h*< hn, 

SO that hn > gn and gn = mm{gn, N'^hn} < g* ■ This implies that n G U^T/. Together, this proves that 
we have properly partitioned [N]. 

We are now ready for the proof of the upper bound on capacity. We argue that the capacity of the 
diamond network with relays cannot be much larger than the sum of the capacities of the L subchannels 
induced by the partition of [N] defined above. Formally, we argue that 

C(iV,(7[^], V]) <C{\T\gT^XLhT^)+C{\T\gT2,2hT^) 

2 L L 

+ + E C{\SuA,gs,,XLhs,,). (23) 

1=1 1=0 k,e=o 

To see this, note that the right-hand side is the capacity of L parallel diamond networks each with unit 
input power constraint. Moreover, increasing each channel gain by a factor of \/2L (or y/2 in the 
case of T^) is equivalent to reducing the power of the additive noise at the destination node of the parallel 
networks by a factor 1/(2L) (or 1/2 for T^). We can now use these parallel networks to simulate the 
original A^-relay diamond network by forcing the input (at the source node u) to all the parallel networks 
to be identical, and by summing up the outputs (at the destination node w) of the parallel networks. This 
proves (l23l) . 

Next, we argue that the capacities of the asymmetric subnetworks in (l23l) can be upper bounded by the 
capacities of symmetric diamond networks. Consider the subset Sk/. Since capacity is increasing in the 
channel gains, 

C{\Sk,il9s,,,,2LhsJ < C(|5fc,,|,2-V,^2i-^/i*). (24) 

Observe that the right-hand side is the capacity of a symmetric diamond network. Consider then T/. By 
the same argument 

C(|T/|,^7Ti,2L/iri) < C(|T/|,2-V,oo), (25) 

and 

C{\T^\,gT2,2LhT2) < C(|T2|,oo,L2i-V). (26) 
It remains to consider and T^. For the set T\ we have 

C{\T'\,gTi,2LhT^) <C{N,N~^g\oo). 

From Theorem [2l 

C{N, iV-V, 00) < 1 log(l + N-'g''). 

By the definition of g*, there exists at least one n such that gn > g* and hn > N~'^g*. Using just this 
one relay t;„, a rate of at least 

ilog(l + iV-V) 

is achievable H For this n, we have 

gn>9''> iv~'^?*, 

^This rate is achievable, for example, with decode-and-forward. Note that we use decode- and-forward here only as a proof technique to 
obtain the upper bound on capacity. Achievability is based exclusively on (bursty) amplify-and-forward. 
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and hence n ^ T^. Moreover, using (|2TI) . 

/in > Ar~V > ^^^^/i'' > N"^h\ 

and hence n ^ T^. This n is therefore an element of one of the subsets {T/}, {T^}, {5*^,^}, and we obtain 
from (EU-dlil), 

C{\T^\,gTi,LhTi) <max|^ max^^C(|T/|, 2-V, oo), ^ max^^ C(|T/|, oo, L2^-^/i*) , 

max C{\S,,,\,2~'g*,L2'~'h'^)]. (27) 

Similarly, 

C{\T'^\,gT2,2hT2) < C{N,oo,2N^^h*) 
< |log(l + 2A^~i/i*) 
<ilog(l + /i*). 

By the definition of h* there exists at least one n such that /i„ > /i* and > h*. Using this relay Vn 
alone, we achieve at least a rate of 

|log(l + /i*). 

For this n, 

hn>h*> N-^h*, 

and hence n ^T'^. Moreover, using (l22l) . 

and hence n ^ T^. This n is therefore an element of one of the subsets {T/}, {T/}, {5*^,^}, and we obtain 
again from (I24l)-(l26l). 

C{\T^\,gT2,2hT2) < maxj^ max^^C(|T/|, 2-V, oo), ^ max^^ C(|T/|, oo, L2^-V), 

max C{\SkA,2~'g\L2'~'h*)]. (28) 
Substituting dSJ-dlH]) into ([23]), we obtain 

C{N,giN],hiN]) <Lmax|^ max^^C(|T/|, 2-V, oo), ^ max^^ C(|T/|, oo, L2^-V) , 

max C{\S,,,\,2~'g*,L2'~'h'^)]. (29) 

This concludes the proof of the upper bound on capacity. 

We continue with the proof of achievability. Fix k,i G {0,...,L}, and recall that a„ is the constant 
determining the amplification at relay Vn- Assume we set a„ = for all n ^ 5*^,^. This results in a network 
in which all but the relays in 5*^,^ are removed. Thus 

R5{N,g[N],h[N]) > Rs{\Sk,i\,gSk,vhs^^,). (30) 

Moreover, since Rs is increasing in the channel gains, 

Rs{\SkA,9s,,,hsJ > Rs{\S,,l2-''-'g\2-'-'h*). (31) 
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With this, we have lower bounded the rate achievable for the asymmetric diamond network by the one 
of a symmetric diamond network (with fewer relays and smaller channel gains). We can thus apply the 
results from Section ITlI-AI to obtain 

sup Rs{\Sk,el2-''-'g\2-'-'h') > ^C(|5fc,,|, 2-V, ^2^-^), (32) 
<5e(o,i] 112L 

where the factor 1/(1 12L) = 1/(8L x 14) is composed of a factor 8L to offset the increase of the channel 
gains to the relay by two and the increase of the channel gains from the relays by AL (see Theorem [T]) 
and of a factor 14 to go from rate achievable with bursty amplify-and-forward to capacity (see Theorem |2] 
and Corollary [3]). 

Combining dlB, and ^ yields 

sup Rs{N,giM],hiN]) > -^C{\SkA,'2-'9\ L2'-'h*). (33) 
<56{o,i] 112L 

A similar argument, setting a„ = for n outside T^, shows that 

sup Rs{N,g[N],h[N]) > sup Rs{\Tl\, gj^i, hrpi) , 
<5e(o,i] 5e(o,i] 

> sup i?4|T/|,2-^-V,2-^-V), 

56(0,1] 

>4c(|T/|,2-V,oo), (34) 



and 



sup Rs{N,g[N],h[N{) > sup i?5(l^/|,5'T2, 

-56(0,1] 56(0,1] 

> sup Rs{\Tl\,N^2~'~^h\2~'~'h*)., 

56(0,1] 



> ic(|T/|,oo,L2i^V), (35) 



56L 

for all £g {0,...,L}. 

We can optimize over the lower bounds in (l33l) . (l34l) . and (l35l) to obtain 



sup Rs(^N, g[N], h[N]j > ^ max | max C(|T/|,2 ^g*,oo), max C(^\T^\,oo,L2^ ^^*), 

56(0,1] 112L lie{0,...,L} £6{0,...,L} 

max C{\SkA,2~'g\L2'~'h*)]. 

fc,£6{0,...,L} ^' ' ' 

Comparing this with the upper bound (|29|) shows that 

C(iV,5([Ar],/l[Ar]) < 112L^ SUp i?5(iV, ^([Ar], /l[Ar]). 

56(0,1] 

Using that 

L < (3 log(A^) + 1)^ + 6 log(iV) + 4 

shows that there exists a universal constant K < oo (and, in particular, independent of gy^^, h^N], and A^) 
such that 112L^ < A'log'^(A^) for N > 2. This concludes the proof of the theorem. ■ 
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V. Conclusion 

We presented an approximation of the capacity of the symmetric Gaussian A^-relay diamond network. 
The capacity was characterized up to a 1.8 bit additive gap and a factor 14 mukiplicative gap uniformly 
for all channel gains and number of relays. The inner bound in this approximate characterization relies 
on bursty amplify-and-forward, showing that this scheme is good simultaneously at low and high rates, 
uniformly in the channel gains and in the number of relays N. The upper bound resulted from a careful 
evaluation of the cut-set bound. We argued that all 2^ possible cuts in the diamond network need to 
be evaluated simultaneously, and that the standard approach of only considering the minimum of the 
broadcast and multiple-access cuts is insufficient to derive uniform capacity approximations. We extended 
this approach to asymmetric diamond networks, for which we showed that bursty amplify-and-forward 
achieves capacity up to a multiplicative gap of a factor 0(log^(A^)) with pre-constant in the order notation 
independent of the channel gains. 

The results in this paper show that, at least for symmetric diamond networks, it is possible to de- 
rive capacity approximations that are independent of the network size. Deriving such uniform capacity 
approximations for general networks remains an open problem. 
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